Noised induced phase transition in an oscillatory system with dynamical traps 
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A new type of noised induced pliase transitions is proposed. It occurs in noisy systems with 
dynamical traps. Dynamical traps are regions in the phase space where the regular "forces" are 
depressed substantially. By way of an example, a simple oscillatory system {x, v = x} with additive 
white noise is considered and its dynamics is analyzed numerically. The dynamical trap region is 
assumed to be located near the x-axis where the "velocity" v of the system becomes sufficiently 
low. The meaning of this assumption is discussed. The observed phase transition is caused by 
the asymmetry in the residence time distribution in the vicinity of zero value "velocity". This 
asymmetry is due to a cooperative effect of the random Langevin "force" in the trap region and the 
regular "force" not changing the direction of action when crossing the trap region. 



I. INTRODUCTION 

The ability of noise to produce order in systems, in 
particular, to induce phase transitions is well established 
(see, e.g., Refs 0,0)- Such a phase transition manifests 
itself in the phase-space density of the system chang- 
ing its structure, for example, the number of maxima. 
Noised-induced phase transitions are distinguished from 
the classical ones by the fact that their cause is not only 
the features of regular "forces" but also the action of 
random Langevin "forces". As a result, in particular, 
the maxima of the distribution function describing the 
noise-induced phases are not necessarily related to the 
zero value of the regular "forces" . 

System of elements with motivated behavior, e.g., fish 
and bird swarms, car ensembles on highways, stock mar- 
kets, etc. often display noise-induced phase transitions. 
The formation of a new phase is caused by noise action 
(for a review see Ref. Q). The theory of these phenom- 
ena is far from being developed well. 

The present Letter considers a certain class of such 
systems whose dynamics can be described by two vari- 
ables X, V which perform a damped harmonic oscillation 
near the equilibrium point {x — 0, v — 0}. However, 
the system "cost" of deviation from the equilibrium can 
differ substantially for these variables. For example, in 
driving a car the control over the relative velocity v is of 
prime importance in comparison with the correction of 
the headway distance x. So, under normal conditions a 
driver should eliminate the relative velocity between her 
car and a car ahead first and only then correct the head- 
way. In markets the deviation from the supply-demand 
equilibrium reflecting in price changes also has to ex- 
hibit faster time variations than, e.g., the production 
cost determined by technology capabilities. In physi- 
cal systems this situation can be also met, e.g., in Pd- 
metal alloys charged with hydrogen where the structure 



relaxation exhibits non-monotonic dynamics P, In 
these alloys hydrogen atoms and nonequilibrium vacan- 
cies form long lived complexes affecting essentially the 
structure relaxation. Their generation and disappear- 
ance governed, in turn, by the structure evolution causes 
the non-monotonic dynamics which can be described in 
terms of dynamical traps. 

These observations lead to the concept of dynamical 
traps, a certain "low" dimensional region in the phase 
space where the main kinetic coefficients specifying the 
characteristic time scales of the system dynamics become 
sufficiently large in comparison with their values outside 
the trap region 0, Q • A trap region is not necessarily 
to be bounded in all the dimensions and, in this case, it 
itself cannot lead to the formation of new phases. More- 
over, the equilibrium point can be absolutely stable, so 
without noise there is only one phase in the system. The 
present Letter analyzes the effect of noise on such a sys- 
tem and demonstrates that additive noise in a system 
with dynamical traps is able to give rise to new phases. 
It should be noted that in most models the phase tran- 
sitions are caused by multiplicative noise, however, addi- 
tive noise also can give rise to these phenomena [sUgHToj. 



II. NOISED OSCILLATORY SYSTEM WITH 
DYNAMICAL TRAPS 

By way of example, the following system typically used 
to describe the harmonic oscillatory dynamics is consid- 
ered 
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FIG. 1: Characteristic structure of tlie pliase space {x,v^. 
Tiie sliadowed domain represents tiie trap region wiiere tlie 
regular "force" is depressed and the system motion is random. 
The regular "force" depression is described by the factor Q.{v) 
illustrated in the left window. The essence of the trap effect on 
the system dynamics is shown in the right window. Outside 
the trap region the system dynamics is mainly regular. 



Here x and v are the dynamical variables usually treated 
as a coordinate and velocity of a certain particle, ujq is the 
circular frequency of oscillations provided the system is 
not affected by other factors, a is the damping decrement, 
and the term eo^u(t) in equation l^l is a random Langevin 
"force" of intensity eo proportional to the white noise 



(3) 



with unit amplitude. The function Q{v) describes the 
dynamical trap effect arising in the vicinity of zero value 
velocity. For this function, the following simple Ansatz 
can be used 
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where the parameter i9t characterizes the thickness of the 
trap region and the parameter A < 1 measures the trap- 
ping efficacy. When A = 1 the dynamical trap effect is 
ignorable, for A = it is most effective. 

The characteristic features of the given system are il- 
lustrated in Fig. n The shadowed domain shows the 
trap region where the regular "force" , the former term in 
Eq. is depressed. The latter is described by the factor 
Q.{v) taking small values in the trap region (for A ^ 1). 
Inside the trap region the system is mainly governed by 
the random Langevin "force" . Outside the trap region it 
is approximately harmonic. 

In order to analyze the system dynamics a dimension- 
less time t and the dynamical variables "q and u are used. 
Namely, the time t is measured in units of I/wqj 
t — !■ t/ijQ and the units of the coordinate x and the ve- 
locity V are 'dt/'^o a-nd i9t, respectively. So, by introducing 



rj = and u 



V 



the dynamical equations (Q), Q read (for the dimension- 
less time t) 

^=u, —^-n[u]{i^ + au) + em. (5) 

where the noise £,{t) obeys conditions like equalities 
the parameter e = to/ {^JujQ'dt), and the function 
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Without noise, this system has only one stationary 
point {t? = 0,M = 0} being stable because it possesses 
a Liapunov function 
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This Liapunov function attains the absolute minimum at 
the point {77 = 0, u = 0} and obeys the inequality 



dH(?7, u) 
Jt 



< for u^Q. 
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In particular, if cr = and e = then function ^ is the 
first integral of the system. In what follows, the values a 
and e will be treated as small parameters. 

The present Letter demonstrates the fact that the noise 
f (t) can cause a phase transition in the given system. It 
manifests itself in that the distribution function V{rj,u) 
changes form from a unimodal to a bimodal one. The dy- 
namics of system ^ was analyzed numerically using a 
high order stochastic Runge-Kutta method ^31 (see also 
Ref. ^3)- The distribution function 7^(77, it) was calcu- 
lated numerically by finding the cumulative time during 
which the system is located inside a given mesh on the 
{rj, u)-plane for a path of a sufficiently long time of mo- 
tion, t K, 500000. The size of mesh was chosen to be 
about 1% of the dimension characterizing the system lo- 
cation on the {rj, it)-plane. 

The evolution of the distribution function Virjju) is 
shown in Fig. |2| in the form of the level contours divid- 
ing the variation scale into ten equal parts. The upper 
window corresponds to the case of A = 1 where the trap 
effect is absent and the distribution function is unimodal. 
The third window illustrates the case when the distri- 
bution function has the well pronounced bimodal shape 
shown also in Fig.|21 Comparing the three upper windows 
in Fig. |2] it becomes evident that there is a certain rela- 
tion $c(A,cr, e) — between the parameters A, cr, and 
e when the system undergoes a second order phase tran- 
sition which manifests itself in the change of the shape 
of the phase space density V(r],u) from unimodal to bi- 
modal. In particular, for cr = 0.1 and e = 0.1 the critical 
value of the parameter A is Ac(cr, e) w 0.5 as seen in the 
second window. 
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FIG. 2: Evolution of the distribution function Virj, u) (shown 
by level contours) as the parameter A decreases. In numerical 
calculations the values a = 0.1 and e = 0.1 were used. The 
lower window depicts only one maximum of the distribution 
function. 
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FIG. 3: The form of the distribution function 'P{ri,u) for the 
parameters a = 0.1, e = 0.1, and A = 0.2. 
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FIG. 4: A typical fragment of the system path going through 
the trap region. The parameters g = 0.1, e = 0.1, and A = 
0.01 were used in numerical simulations in order to make the 
trap effect more pronounced. 



III. MECHANISM OF THE PHASE 
TRANSITION 

To understand the mechanism of the noised induced 
phase transition observed numerically in the given sys- 
tem, consider a typical fragment of the system motion 
through the trap region for A <C 1 that is shown in Fig.^] 
When it goes into the trap region Qt, —'dt ^ w <C "i?*, the 
regular "force" ^[u] {rj + au) containing the trap factor 
f2[u] and governing the regular motion becomes small. So 
inside this region the system dynamics becomes random 
due to the remaining weak Langevin "force" {t) . How- 
ever, the boundaries d^Qt (where v ~ 'dt) and d-Qt 
(where v ~ — ^9t) are not identical in properties with 
respect to the system motion. At the boundary d+Qt 
the regular "force" leads the system inwards the trap re- 
gion Qt, whereas at the boundary d-Qt it causes the 
system to leave the region Qt- Outside the trap region 
Qt the regular "force" is dominant. Thereby, from the 
standpoint of the system motion inside the region Qt, 
the boundary d+ Qt is "reflecting" whereas the boundary 
d-Qt is "absorbing". 

As a result the distribution of the residence time at 
different points of the region Qt should be asymmet- 
ric, as schematically shown in Fig. ^the right window). 
This asymmetry is also seen in the distribution function 
V{ri,u) obtained numerically. Its maxima are located at 
the points with non-zero values of the velocity, which is 
clearly visible in the lower window of Fig. |21 Therefore, 
during location inside the trap region the mean velocity 
of the system must be positive and it tends to go away 
from the origin. This effect gives rise to an increase in the 
"energy" 1-L{rj,u). Outside the trap region the "energy" 
T-l{r],u) decreases according to expression Q. So, when 
the former effect becomes sufficiently strong, i.e., the ran- 
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dom "force" intensity e exceeds a certain critical value, 
e > ec(A,cr), the distribution function V{'q,u) becomes 
bimodal. 

The system location with respect to the velocity v is 
due to the regular "force" being sufficiently strong out- 
side the trap region, so the system spends the main time 
inside this region. Its location with respect to the co- 
ordinate X is caused by the fact that the region where 
the Langevin "force" mainly affects the system dynamics 
decreases in thickness as the coordinate x increases. The 
latter tendency takes place because the regular "force", 
the first term in Eq. is proportional to x. In fact, 
the thickness U{7f} of the trap region in the vicinity of the 
point {77, u — 0\ can be estimated using the condition of 
the equality of the characteristic times, tg and td, during 
which the system crosses the trap region under action of 
the regular "force" and the random Langevin "force" . So 

^ r![[/]77 ~ ~ £2 ■ 

and setting for the sake of simplicity A = we get the 
estimate 

C/(?7) ~ e'/^-i/^ 

Moreover, let the mean velocity in the trap region caused 
by the residence time asymmetry be about U{ri). Then 
the characteristic increase Srj of the coordinate 77 got by 
the system when crossing the trap region is estimated 
as Srj ^ Utg l/rj. Thereby the dynamical trap effect 
becomes weaker as the "energy" 7i(??, u) increases. By 
contrast, according to Exp. iQ the higher the "energy", 
the stronger its dissipation caused by the regular "force" . 



IV. CONCLUSION 

The present paper has considered noise-induced phase 
transitions in systems with dynamical traps. By way 
of example, a simple oscillatory system {x, v = x} is 
studied when the trap region is located in the vicinity of 
the X-axis and without noise the stationary point {x = 
0,v ~ 0} is absolutely stable. For this system as shown 
numerically the additive white noise can cause the phase- 
space density to take the bimodal shape. 

In contrast to the classical phase transitions the posi- 
tion of the new noise-induced phases is not specified by 
the zero- values of regular "forces" even approximately. 
The cause of the observed phase transition is the asym- 
metry of the residence time distribution inside the trap 
region. This asymmetry is due the cooperative effect of 
the regular "force" outside the trap region and the ran- 
dom Langevin "force" inside it. The regular "force" does 
not change the direction when crossing the trap region, 
inside this region it is depressed only. As a result, for 
the motion inside the trap region one of its boundaries is 
"reflecting" , whereas the other is "absorbing" , which in- 
duces the residence time asymmetry. The latter gives rise 
to increase in the system "energy" . Outside the trap re- 
gion the regular "force" causes the "energy" to decrease. 
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